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Abstract We propose and theoretically investigate an 
unambiguous Bell measurement of atomic qubits as¬ 
sisted by multiphoton states. The atoms interact res¬ 
onantly with the electromagnetic field inside two spa¬ 
tially separated optical cavities in a Ramsey-type in¬ 
teraction sequence. The qubit states are postselected 
by measuring the photonic states inside the resonators. 
We show that if one is able to project the photonic 
field onto two coherent states on opposite sites of phase 
space, an unambiguous Bell measurement can be im¬ 
plemented. Thus our proposal may provide a core el¬ 
ement for future components of quantum information 
technology such as a quantum repeater based on coher¬ 
ent multiphoton states, atomic qubits and matter-field 
interaction. 


1 Introduction 

Establishing well-controlled entanglement between spa¬ 
tially separated quantum systems is essential for quan¬ 
tum communication m- At its core a quantum re¬ 
peater employs entanglement which is generated and 
distributed among intermediary nodes positioned not 
too distant from each other. Entanglement purifica¬ 
tion m enables the distillation of a high-fidelity state 
from a large number of low-fidelity entangled pairs and 
with the help of entanglement swapping procedures 
the two end points of a repeater are entangled. There 
are many different implementation proposals for quan¬ 
tum repeaters, utilizing completely different systems 
and entanglement distribution protocols [5]. A promis- 
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ing approach towards these schemes is to require some 
compatibility with existing optical communication net¬ 
works. The proposal of van Loock et al. [7l[8l[^IT(1] is 
such an approach where the repeater scheme employs 
coherent multiphoton states. These proposals assume 
dispersive interaction between the atomic qubits and 
the single-mode of the radiation field. This imposes 
limitations on the photonic postselection. It was shown 
that these limitations can be overcome in the case of res¬ 
onant atom-field interactions mm and it was demon¬ 
strated for one building block of a repeater, namely the 
entanglement generation between spatially separated 
and neighbouring nodes. A natural extension of this 
approach is to propose resonant atom-field interaction 
based schemes also for the other building blocks. In the 
case of entanglement swapping a complete atomic Bell 
measurement is required. 

Bell measurements play a central role also in en¬ 
tanglement-assisted quantum teleportation m and in 
superdense coding [TS]. In the case of photonic qubits 
theoretical proposals miiiziiiH] have been made and ex¬ 
perimental realizations have already been carried out 
[Hun]. However, for atomic qubits there are still ex¬ 
perimental difficulties which restrain implementations 
of complete Bell measurements where projections onto 
the four Bell states can be accomplished. There exist 
experimental proposals that rely on the application of a 
controlled NOT gate These proposals have the 

drawback that experimental implementations of two- 
qubit gates have still complications to attain high fi¬ 
delity This implies that the fidelity of the 

generated Bell states is also affected [53]. A proposal 
focusing specifically on a non-invasive atomic Bell mea¬ 
surement with high fidelity is still missing. 

In previous work we have introduced a protocol to 
project onto one Bell state with high fidelity [T3] based 
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on atomic qubits which interact sequentially with co¬ 
herent field states prepared in two cavities. The field 
states emerging after the interactions are postselected 
by balanced homodyne photodetection. In this paper 
we expand our previous work to accomplish the pro¬ 
jection onto all four Bell states provided the protocol 
is successful. Thus we introduce an unambiguous Bell 
measurement of two atomic qubits with the help of co¬ 
herent multiphoton field states. We demonstrate that 
the possibility of implementing field projections onto 
two coherent states on opposite sites of phase space 
implies the possibility to realize an unambiguous Bell 
measurement. Our protocol has a finite probability of 
error depending on the initial states of the atoms. This 
is due to the imperfect overlap of the field contributions 
with coherent states. Nevertheless, it is an unambiguous 
protocol as there are four successful events that lead to 
postselection of four different Bell states. The scheme 
is based on basic properties of the two-atom Tavis- 
Cummings model |26j and on resonant matter-field in¬ 
teractions which are already under experimental inves¬ 
tigation These considerations make our 

scheme compatible with a quantum repeater or a quan¬ 
tum relay based on coherent multiphoton states, atomic 
qubits and resonant matter-field interaction. Our pro¬ 
posal demonstrates that scenarios involving the two- 
atom Tavis-Cummings model are rich enough to enable 
future Bell measurement implementations. 

The paper is organized as follows. In Sec. we in¬ 
troduce the theoretical model and analyse the solutions 
of the field state with the aid of the Wigner function 
in phase space. Furthermore, we provide approximate 
solutions of the global time dependent state vector that 
facilitate the analysis of the system. In Sec.j^we present 
a scheme to perform an unambiguous Bell measurement 
provided one is able to project a single mode photonic 
field onto coherent states. In Sec. we provide a numer¬ 
ical analysis of the fidelity of the projected Bell states 
and discuss general features of the protocol. Details of 
our calculations are presented in Appendices and 


2 Theoretical model 

2.1 Basic equations 

In this section we recapitulate basic features of the 
two-atom Tavis-Cummings model [26]. This model has 
been considered previously to study the dynamics of en¬ 
tanglement [OII551I51II351I55] . The model describes the 
interaction between two atoms A and B and a sin¬ 
gle mode of the radiation field with frequency oj. The 
two identical atoms have ground states |0)i and excited 


states |I)i (i € {A,B}) separated by an energy differ¬ 
ence of Huj. In the dipole and rotating-wave approxima¬ 
tion the Hamiltonian in the interaction picture is given 

by 

H = hg {ata + a-a'f) (1) 

i=A,B 

where o-j*" = |I) (0|i and d” = |0)(I|i are the atomic rais¬ 
ing and lowering operators {i G {A,B}), and a (a’^) is 
the annihilation (creation) operator of the single mode 
field. The coupling between the atoms and the field is 
characterized by the vacuum Rabi frequency 2g. 

The time evolution of the system can be evaluated 
for an initial pure state as 

( 2 ) 

We are interested in the case where the atoms and the 
cavity are assumed to be prepared in the product state 

\rp,) = (^c-\^-) + c+\^+)+d^\^^)+d+\^+))\a), (3) 

with the radiation field considered initially in a coherent 
state 

2 

|a) = ^e-^4=|n), a = (4) 

with mean photon number n and photon-number states 
|n). The parameters and are the initial probabil¬ 
ity amplitudes of the orthonormal Bell states 

= (5) 

with the atomic states \i,j) = \i)A\j)B G {0,1}). 
We have chosen an atomic orthonormal basis containing 
the states as the state \'I'~)\n) is an invariant state 
of the system. This is explained in Appendix [^ where 
we present the full solution of the temporal state vec¬ 
tor. The other two Bell states \d>^) depend on the initial 
phase of the coherent state. They appear naturally 
in the Tavis-Cummings model due to the exchange of 
excitations between atoms and cavity, and are involved 
in an approximate solution of the state vector that facil¬ 
itate the analysis of the dynamics. Before showing the 
detailed form of our solution, let us give an overview of 
the dynamical features that impose relevant time scales 
in the system. 
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2.2 Collapse and revival phenomena 

The collapse and revival phenomena of the Jaynes-Cum- 
mings model and of the two-atom Tavis-Cummings model 
play an essential role in the quantum information proto¬ 
cols presented in Refs. nniniEa. This behaviour was 
first found in the time dependent atomic population 
in the Jaynes-Cummings model m when the field is 
initially prepared in a coherent state: The populations 
display Rabi oscillations that cease after a collapse time 
tc and appear again at a revival time tr- In the case of 
the two atom Tavis-Cummings model, the collapse and 
revival time of the Rabi oscillations are given by 

tj. = —^jAn + 2, tc = 

9 

These time scales have been previously introduced and 
can be found, for instance, in Refs. |33l|37]. As they 
play an essential role in the dynamics of the system it 
is convenient to introduce the rescaled time 

T = t/tj. = tg/iT^/ An + 2. (7) 

Let us explain these phenomena by visualizing the 
phase space of the radiation field with the aid of the 
Wigner function [38ll^ 

= ^J ( 8 ) 

with the complex numbers /3 and (^. The operator gt = 
Tr atoms {l^t) (^t 1} is the density matrix of the field state 
obtained after taking partial trace over the atomic de¬ 
grees of freedom from the full density matrix corre¬ 
sponding to the state vector in Eq. (§. In Fig. we 
show the Wigner function after interaction times r = 
1/4 in the top panel and r = 1/2 in the bottom panel. 
The circular shape corresponds to the initial coherent 
state I a). This contribution to the field remains station¬ 
ary as long as there is an initial contribution of the state 
The reason is that \’I'~)\n) is an invariant state 
of the system. There are two other contributions to the 
field that rotate around the origin. In the top panel of 
Fig. 0 it can be noticed that for an interaction time 
of r/4 they have completed a quarter of cycle. At the 
bottom, the situation at interaction time r/2 is shown 
where half a rotation has been completed. The interfer¬ 
ence fringes between the field contributions signify that 
there are coherent superposition between these states 
of the field. The behaviour of the field state in phase 
space explain the phenomena: Rabi oscillations cease 
(collapse) when the field contributions are well sepa¬ 
rated, e.g. at time t/4, and revive when the field con¬ 
tributions overlap, e.g. at r/2 or the main revival at r 
when all the field constituents coincide at the position 
of the initial coherent state. 
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Fig. 1 Wigner function of the cavity field after the inter¬ 
action with two two-level atoms at times r = 1/4 (top) and 
T = 1/2 (bottom) with the rescaled time of Eq. 0 . The 
initial states of the atoms are defined by the parameters: 
c- = 0.5554, c+ = 0.3213 -|- *0.5004, d.- = -0.2053 -|- *0.3726, 
dj = 0.1046 -|- *0.3819, and the parameter a = \/36.16e*^'®^ 
characterizes the initial coherent sate. 


2.3 Approximation of the state vector 


The full solution to the time dependent state vector 
of the two atoms Tavis-Cummings model has already 
been presented in previous work, see for instance |351 
HD]. Coherent state approximations have also been con¬ 
sidered in the past [I3IM1IM1I1I]. In this context, the 
eigenfrequencies of the Hamiltonian Q that depend on 
the photonic number n are expanded in a first order 
Taylor series around the mean photon number n. How¬ 
ever, the coherent state description is accurate only for 
times well below the revival time. In this work we go 
beyond the coherent state approximation by consider¬ 
ing second order contributions of the eigenfrequencies 
around n. The details can be found in the Appendix [A| 
where it is shown that the time dependent state vector 
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of the system can be approximated by 

|a) + 

(|‘f'+) - l'?J+2..)) I«x) + 
(1^^+) + |<?J_2.r)) l«;), 

with the photonic states 


= E^ 

n=0 


±i 27 rr n+l+n— 




and with the normalization factor 



Fig. 2 Fidelity of the total state of Eq. § with respect to 
the exact solution given by Eq. ( |20[ | as a function of the time 
T in Eq. 0 scaled in terms of the revival time tr'. Five curves 
are presented for different values of the mean photon number 
n as described in the legend. The rest of the parameters are 
the same as in Fig. ^ 


= (l + Re[{c+ + dj)*(c+ - |a+)] sin^(27rT) 

+ 2Ile[d^(dj)*]lm[(a- |a)] sin(27rT). 

\ 1/2 

+ 2Im[(c+)*d^]Re[(a“|a)] sin(27rr)j . (11) 


The quantity {a^ |a) is evaluated in Appendix [b] and 
an approximate expression is given in Eq. (29). 

In order to test the validity of Eq. ®> we have con¬ 
sidered the fidelity F{t) = \{F^\Ft^r)\‘^ of the approxi¬ 
mated state vector with respect to the exact result given 
in Eq. (20). In Fig. we have plotted the results of 
numerical evaluations of the fidelity F{t) for different 
values of the mean photon number h. It can be noticed 
that the validity of this approximation improves with 
increasing mean photon number n. In the Appendix [A| 
it is discussed that our approximation is valid provided 
the condition r ^ ^/nj2TT is fulfilled. 

The form of the solution given in Eq. (§ allows a 
simple analysis of the dynamics. It is written in terms 
of an orthonormal atomic basis of Bell states and is 
therefore suitable for the analysis of the atomic entan¬ 
glement. In particular, it is interesting to note that for 
an initial state without a contribution of the state |^^), 
i.e. = 0, a photonic projection that discriminates the 
state I a) from the states |a^) can postselect the atomic 
Bell state \F~). In Ref. [13] we studied this Bell state 
projection and found that its implementation requires a 
flexible restriction for the interaction time: it has to be 
below the revival time and above the collapse time given 
in Eq. §)• In the following we concentrate in a more 
specific interaction time. We analyze the dynamics at 
the specific interaction time r = 1/2. This analysis will 
allow us to introduce in Sec. a protocol to perform 
the four Bell state projections. 


2.4 Basic dynamical features at scaled interaction time 
r = 1/2 

There are two main reasons for studying in detail the 
case with scaled interaction time t = 1/2. The first 
one is that the time dependent atomic states in Eq. ([^ 
coincide, i.e. 


J = -\^p. (12) 

The second reason is that the photonic states |Q!j*/ 2 ) 
have completed half a rotation in phase space and lie 
on the opposite site to the initial coherent state |a) 
whereby overlapping with the coherent state |—o). This 
means that at this interaction time and for |a| 3> 1, the 
initial photonic state |a) can be approximately distin¬ 
guished from the other two states However, the 

states overlap significantly. This can be noticed 

in Fig. where we have plotted the Wigner function. 
The circular shape corresponds to the initial coherent 
state la), while the distorted ellipses on the opposite 
site of the phase space correspond to the states 
To distinguish these two components of the field it is 
convenient to conceive an experiment that is able to 
project the field state onto the coherent states |±a). 

In order to study the projection onto the state |±a), 
one has to evaluate its overlap with the photonic states 
of the state vector in Eq. (§. First we consider the 
overlaps that can be neglected for large value h, namely 


l-a) = 


± \ 


, \a\a-^i 2 j C)C e 


(13) 


The explicit form of the overlap {a\a^i^ is given in 
Eq. (28) of the Appendix [b] where its approximation is 
also evaluated. The nonvanishing overlaps in the limit 
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of large mean photon number are (Qf|a) = 1 and 




Vd + TT- 




The expression in Eq. (14) is also evaluated in detail 


in Appendix]^ This overlap is real valued if the mean 
photon number n = |ap fulfills the relation 


n = m + arctan with m S N. 


(15) 


If the condition in Eq. (151 is fulfilled and if we sup¬ 


pose an initial atomic state with no contribution from 
the state i.e. dj = 0, then it can be verified that 

a projection onto the field state |—a) postselects the 
atoms in the unnormalized atomic Bell state \/bc^\^'^), 
with 


6 = 2 / \/4 + TT^. 


(16) 


The success probability of this projection is given by 
6|c+p which is proportional to the initial probability of 
this particular Bell state The factor 6 is the result 
of our inability to project perfectly and simultaneously 
onto both field states section we 

present a protocol that can perform postselection of 
the four Bell states regardless of the initial state of the 
atoms. 


3 An unambiguous Bell measurement 

In this section we introduce a protocol which imple¬ 
ments a projection onto four orthogonal atomic Bell 
states of Eq. for any given initial condition of the 
atoms. The scheme we propose requires interactions be¬ 
tween the atoms with two different cavities as sketched 
in Fig. The interaction time between the atoms and 
the electromagnetic field in each cavity is assumed to 
be T = 1/2. The field in the first (second) cavity has 
to be prepared in a coherent state |a) (|*a)). After the 
interaction with the first cavity the resulting field is 
projected onto the initial state |a). In case of failure a 
projection onto the state |—a) is performed. The pro¬ 
jection of the field postselects the atoms in a state that 
has contribution of only two of the Bell states. This 
postselected atomic state is taken as initial condition 
to interact with a second cavity prepared in the state 
\ia). The atoms are assumed to evolve freely for a time 
Tf before interacting with a second cavity. This does not 
affect the protocol as the free Hamiltonian commutes 
with the interaction Hamiltonian in Eq. 0 . After the 
interaction of the atoms with the second cavity, the 
field in the second cavity is projected onto \ia) and if 
this fails another projection onto the state |—fa) is per¬ 
formed. With this field state projection, the atoms are 



I . T. I'l .T f . t. T .I 


Fig. 3 Schematic representation of the proposed atomic Bell 
measurement: Two atomic qubits interact with the electro¬ 
magnetic field inside two independent cavities in a Ramsey- 
type interaction sequence. Different projections on the field 
states inside the cavities, recorded by detectors Di and D 2 , 
result in a postselection of atomic Bell states as described in 
Table □ 

finally postselected in a unit fidelity Bell state. In what 
follows we discuss in detail all the possible outcomes of 
the protocol. There is a finite probability to fail com¬ 
pletely when none of the coherent state field projections 
is successful. This is discussed in Sec. IH 


3.1 Projection onto |a) in the first cavity 


Let us consider a successful projection onto the field 
state la) of the first cavity. In this case the atoms are 
postselected in the state 

I 

7K 




(17) 


with probability Pi = |c \‘^ + \d^\'^. To write this state 
we have also considered the relations 


1^.- 


<p+Tr/ 2 ') * 1 ^ 0)1 ^^+Tr /2 


= idi. 


(18) 


The postselected atomic state of Eq. 0 is taken 
as initial condition for the interaction with the second 
cavity prepared in the coherent state \ia) as depicted 
in Fig.[^ Two scenarios are possible for the projection 
of the field in the second cavity. In the first place we 
consider a projection onto the coherent state \ia) where 
the atoms are postselected in the state with prob¬ 
ability Pii = |c“p/Pi. This can be verified from Eq. 
0 as the new initial state does not have a contribu¬ 
tion of projections performed in the 

first and second cavity are independent events the state 
\d'~) can be projected with overall success probabil¬ 
ity PiPii = |c“p, the initial probability weight of this 
state before the protocol. The second possibility is to 
project onto the state \—ia). In that case the atoms are 
postselected in the state provided 

the condition in Eq. (15) is fulfilled. This can be veri¬ 


fied using Eq. 0 with an initial coherent state \ia) and 
the atoms initially in the state of Eq. p7|) that has no 
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Field state 
in detector 

Di 

Field state 
in detector 
D 2 

Atomic 

state 

|Bell> 

Probability 

|a> 

|iQ> 


|C-|2 

|a> 

\-ia) 



|-o> 

\ia) 



|-o> 

\-ia) 

\F+) 

b2|c+|2 


Table 1 Summary of the Bell state protocol assisted by pho¬ 
tonic state measurements. The first (second) column indicates 
the photonic field that has to be selected in the first (second) 
cavity by detector Di (D 2 ) in the interaction sequence de¬ 
picted in Fig. The third column indicates the resulting 
atomic state with the probability of occurrence given in the 
last column with b = 2/^4 -|- ;:s! 0.537. The protocol fails 

with probability (1 - b){\d-\^ + \d+\^) + {1 - b^}\c+\\ 



n 

Fig. 4 Fidelity Fb of the projected atomic Bell states as a 
function of the initial mean photon number of the fields inside 
the cavities: The interaction time in both cavities is given by 
T = 1/2, i.e. half the revival time (see Eq. 0 ) and the rest 
of the initial conditions are the same as in Fig.[^ Each curve 
correspond to a different Bell state as explained in the legend. 


contribution of I'?"'’). The success probability for this 
event is Pio = b\d'^\‘^/Pi. Correspondingly the projec¬ 
tion onto the atomic state occurs with overall suc¬ 
cess probability PiPio = b\d'^\‘^. This is proportional to 
its initial probability weight but not equal. The propor¬ 


tionality factor b is given in Eq. (16) and accounts to 


the imperfect projection onto the states 

3.2 Projection onto |—a) in the first cavity 

Now we consider a successful projection onto the coher¬ 
ent state |—a) in the first cavity. In this situation the 
atoms are postselected in the state 


^{c+\W+)-d 


0+7r/2 IV 7r/2 


(19) 


with probability Pq = + b\d'\\’^. We have used the 


relations in Eqs. (12) and 


The normalized state of Eq. (19) is taken as initial 


condition to interact with the second cavity prepared in 
the coherent state \ia). There are two scenarios in the 
projection of the second cavity. First we consider a suc¬ 
cessful projection onto the state \ia). As the initial state 


of Eq. (191 does not have any contribution of |<f' ), the 
atoms are postselected in the state 
This occurs with success probability Poi = b\d+\^/P^. 
Thus the state is postselected with an overall suc¬ 
cess probability PqPqi = b\d'^\^. A second possible situ¬ 
ation is a projection onto the state \—ioL) in the second 
cavity. In this situation the atoms are postselected in 
the state 1^'''"). This can be noted from Eq. ([^ as the 
second initial atomic state of Eq. (19) does not have any 


contribution of the state 0 +^/ 2 )- success proba¬ 
bility of this event is Pqo = |c-i-pI'^/Po- It implies an 
overall success probability of postselecting state 


4 Discussion of the protocol 

4.1 Fidelity of the postselected Bell states 

In order to test our protocol based on the approxima¬ 
tions of Eqs. (|^ and (14) we have numerically evalu¬ 
ated the fidelity Pb = | (Belllj//) p of the resulting Bell 
states in each of the four possible successful outcomes. 
The state |Bell) stands for any of the four Bell states 
of Eq. (§. The state |^/>) is the exact numerical solu¬ 
tion after the protocol and depends either on the n or 
r. In Fig. we have plotted the fidelity Pb for the 
different Bell states as a function of the mean pho¬ 
ton number fi = jap of the initial coherent field states 
I a) and \ia). Interestingly, the protocol already shows 
high fidelity (above 0.9) even for small mean photon 
numbers. The results improve for increasing values of 
h in accordance to the validity of our approximation 
for high photon number explained in Sec. The fi¬ 
delity has an oscillatory periodic behaviour and max¬ 
ima are achieved close to the values of h predicted by 


ofPoPoo = 62|c+|2. 


Eq. (15), i.e. when h is an integer number plus the con¬ 
stant arctan(7r/2)/27r « 0.16. A possible error 5n in the 
previous value has to fulfill the condition Sn I/tt to 
ensure a high fidelity of the atomic states. It should be 
mentioned that in the extreme opposite case in which 
both cavities are initially prepared in the vacuum state, 
i.e. fi = 0, the proposed protocol does not work. Ac¬ 
cording to Eq. § the four orthogonal Bell states are 
paired up with three field states and in order to filter 
out all Bell states they have to be orthogonal. This re¬ 
quirement can only be fulfilled in the limit of high mean 
number of photons. 

To test the sensitivity of the protocol with respect 
to the interaction time we have also evaluated the fi¬ 
delity Pb as a function of the scaled interaction time 
T between the atoms and the cavities. The results with 
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Fig. 5 Fidelity Fb of the projected atomic Bell states as a 
function of the scaled interaction time r (see Eq. 0 ) on both 
cavities: The initial conditions are the same as in Fig. ^ with 
n = 36.16. Each curve correspond to a different Bell state of 
Eq.|5] as explained in the legend. 


the initial atomic conditions of Fig. are plotted in Fig. 

We present the results for an initial coherent state 
with mean photon number h = 36 + arctan(7r/2)/27r. 
The black solid curve represents the fidelity of project¬ 
ing onto state \^~) and it shows a constant unit fidelity 
in the time interval of the plot. The stability of this re¬ 
sult has also been discussed in Ref. m and is due to 
the fact that is a special invariant atomic state 
of the two-atom Tavis-Cummings model. The fidelity 
of the state also shows robustness with respect to 
the interaction time r. This is due to the fact that this 
state is obtained after projecting onto \ia) which is the 
stationary initial state of the second cavity. Advantages 
of the two-atom Tavis-Cummings model for generating 
this particular Bell state have also been mentioned pre¬ 
viously in Ref. [Mj. The other two fidelities of projec¬ 
tions onto states and |<F+) oscillate as a function of 
r. In this case the second field projection is performed 
onto field state \ —ia) and this in turn has to “catch” the 
time dependent states \ia^)- Therefore, the oscillations 
are originated by the overlap between photonic states 
{—a\af) that is calculated in the Appendix]^ One can 
estimate that the fidelity Fp around r = f /2 oscillates 
with frequency 2(h -|- f). The optimal interaction time 
according to Eq. (291 is r = 1/2 where the absolute 
value of the overlap attains its maximum. A possible 
error e in the scaled interaction time r = 1/2 -|- e has 
to be restricted to the condition jej ^ l/ 47 r(h -|- 1). 


4.2 Experimental constraints 

Our protocol requires that the pair of atoms interact 
with two different coherent states. This could be real¬ 
ized, for instance, by transporting and positioning the 
atomic qubits in separate cavities. Current experimen¬ 
tal realizations report coherent transport and controlled 
positioning of neutral atoms in optical cavities mm 


mM: where a dipole trap is used as a conveyor belt 
to displace them. Two trapped ions have also been re¬ 
ported to be coupled in a controlled way to an optical 
resonator mM- The cavity can be shifted with respect 
to ions, allowing to tune the coupling strength between 
ions and optical cavity. In this setting, instead of trans¬ 
porting the atoms to a different cavity, the same cavity 
might be shifted to a position where it decouples from 
the atoms until the measurement is achieved. Then it 
would have to be prepared and shifted again for a sec¬ 
ond interaction with the ions. 

In our discussion we have not considered losses. The 
effects of decoherence can be neglected in the strong 
coupling regime, where the coupling strength g between 
atoms and cavity is much larger than the spontaneous 
decay rate of the atoms 7 and the photon decay rate of 
the cavity k. Actually, in our setting due to the specific 
interaction time tr/2 « TT^/Ejg tighter constraints are 
required. More specifically, for the cavities we require 
I/k 3> 'K\Jfilg and for the atoms 1 / 7 ^ ir^/^jg. The 
experiment by Khudaverdyan et al. |42j achieved ratios 
g/n = 32.5 and g /7 = 5 which imply that h ^ 2.5. 
For a single atom interacting with a cavity, the experi¬ 
ment by Birnbaum et al. |45j involves ratios g/n = 8.26, 
gjj — 13.03 and if there is a possibility to attain these 
parameters for a two-atom scenario then the constraint 
would yield h <?; 7. In microwave cavities [44] , the num¬ 
bers are g/K « 60 and 5/7 « 3000 which lead to the 
condition h ^ 360. Thus, the coherence requirement 
of our proposal is in the reach of current experimental 
capabilities. 

We have mentioned that our protocol requires the 
implementation of projections onto coherent states. We 
are not aware of an experimental solution to this prob¬ 
lem. However, coherent states and the vacuum state 
are routinely distinguished in current experiments, see 
e.g. [46]. A successful measurement of the vacuum state 
is achieved when no photons are detected. Therefore, 
for our purposes it would be sufficient to displace the 
state of the field in such a way that the field contribu¬ 
tions |q!^/ 2 ) close to the vacuum state. This can be 
achieved by driving the optical cavity with a resonant 
laser. The Hamiltonian describing this situation in the 
interaction picture is H = h{f2*a + f2a^). Under this 
interaction the states of the field evolve under the in¬ 
fluence of the evolution operator Ut^ = exp {—itdV/h) 
that can be identified with the displacement opera¬ 
tor D{a.) = exp {aa^ — a*a) provided the interaction 
strength of the laser 12 and the driving time are ad¬ 
justed as fitd = ia. In this way one is able perform the 
displacement D{a)\—a) = |0). Finally, we conceive a 
photodetection of the field with three possible outputs; 
I) No signal, meaning a projection onto the vacuum 
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state, i.e. |—a) in the undisplaced picture; 2) a weak 
signal indicating a failure of the protocol; 3) a strong 
signal would come from the field state |2 q!). 


4.3 Probabilities in the protocol 

A summary of all the possible outcomes of the proto¬ 
col is given in the Table We note that summing the 
probabilities of all the successful outcomes of the pro¬ 
tocol results in an overall success probability of Pt = 
6 -I- (1 — 5)(|c“p — 6|c+p) which depends on the ini¬ 
tial state of the system. The complementary probabil¬ 
ity 1 — Pt corresponds to events that lead to failure 
of the protocol. There is a possible failure after a suc¬ 
cessful projection onto ja) but unsuccessful projection 
onto \—ia). This occurs with probability (1 — h)\d~^\'^. 
It also might happen that the projection onto the field 
state |—a) in the hrst cavity is unsuccessful. This takes 
place with probability (1 — 6)(|c+p -I- IdJP). Finally, 
it is possible that both projections in the first and 
second cavity fail with probability (1 — b)b\c'^\'^. Sum¬ 
ming all these failure probabilities leads to 1 — Pt = 
(l-6)(|d-p + M+p) + (l-52)|c+p. 


5 Conclusion 

We have presented a proposal of an unambiguous Bell 
measurement on two atomic qubits with almost unit- 
fidelity. The theoretical description of the scheme in¬ 
volves the resonant two-atom Tavis-Cummings model 
and a Ramsey-type sequential interaction of both atoms 
with single modes of the electromagnetic field in two 
spatially separated cavities. The first and second cav¬ 
ities are initially prepared in coherent states |a) and 
I fa) respectively. The interaction time can be adjusted 
by controlling the velocities of the two atoms passing 
trough the cavities. Our discussion has concentrated 
on basic properties of the two-atom Tavis-Cummings 
model in the limit of high photon numbers. We have 
derived an approximate solution of the dynamical equa¬ 
tion which is expressed as a sum of three terms corre¬ 
lating atomic and field states. A superposition of two 
atomic Bell states is correlated with the initial coher¬ 
ent state. Superpositions of the other two Bell states 
are correlated with two time dependent field states. In 
phase space these time dependent contributions of the 
field state overlap on the opposite site to the initial co¬ 
herent state I a) (|fa)) in the first (second) cavity at an 
interaction time of half the revival time. For this rea¬ 
son we have proposed projections onto the two coher¬ 
ent states |a) and |—a) in the hrst cavity, and \ia) and 
\—ioi) in the second cavity. In order to obtain almost 


unit hdelity atomic Bell states the mean photon num¬ 
ber has to be restricted to the condition given in Eq. 
(15). Our protocol has a hnite error probability due to 
the imperfect projection onto the time dependent con¬ 
tributions of the held states in the cavities that overlap 
with |—a) and \—ia). Nevertheless, the four successful 
events of our protocol summarized in Table unam¬ 
biguously project onto four different Bell states with 
almost unit hdelity. 

In view of current experimental realizations of quan¬ 
tum information protocols in the held of cavity quan¬ 
tum electrodynamics the scheme discussed in this work 
requires cutting edge technology. An experimental im¬ 
plementation would require accurate control of the in¬ 
teraction time and of the average number of photons 
in the cavity. Furthermore, the coherent evolution of 
the joint system must be preserved. This imposes the 
condition that the characteristic time of photon damp¬ 
ing in the cavity and of atomic decay have to be much 
larger than the interaction time that scales with the 
square root of the mean photon number in the cav¬ 
ity. Finally, we point out that the implementation of a 
von Neumann coherent state projection is, up to our 
knowledge, an open problem that has to be considered 
in future investigations. If these obstacles are overcome, 
our proposal offers a key component for quantum infor¬ 
mation technology such as a multiphoton based hybrid 
quantum repeater. 
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A Approximations with large mean photon 
numbers 


In this Appendix we present the derivation of the time depen¬ 
dent state vector of Eq. @. It has been shown in Ref. ESI 
1401 that the time evolution of any initial state in the form of 
Eq. (§ can be obtained from the solution of the eigenvalue 
problem of the two-atom Tavis-Cummings Hamiltonian 0- 
The exact solution can be written in the following form 

m = |0, 0>|x?> + 11, l>|x(> + l^+>IXt+> + c-|<?'->|a> (20) 

with the relevant photonic states 


lx?> = copo|0> -b ^ 




\/2n-l 


h>. 


Ix(> = 




^ 271-1 


|n-2>, 


OO 

IXr+> = E 


(21) 
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and with the aid of the following abbreviations 


= - 

^n.t 




[ ^+Pn—1 T 


y/n CQPn + \/n — l cip„ 


V2n-i 


= - 


y/n-l CoPn - \/ncip „_2 


V2n - 1 


= gy/4n — 2. 


+ 

^3 


02 



The coefficients pn are initial probability amplitudes of the 
photon number states |n) of the initial field state \a). The 
coefficients cq and ci are the initial probability amplitudes of 
the states |0, 0) and |1,1) and are related to the probability 
amplitudes of the state in Eq. ^ by 


Coe*''' ± Cie 


( 22 ) 


The expressions of Eq. ( |21[ | can be significantly simpli¬ 
fied approximately by taking into account that the field is 
initially prepared in a coherent state |a} with photonic distri¬ 
bution p„ = exp(—h/2 -|- i(j>) n'^ /n\ and by assuming a large 

mean photon number n = |ap ^ 1. In such case the photonic 
distribution has the following property 


Pn = \l ~ e^^pn-l. 

n 


(23) 


Applying this approximation to the states of Eq. ( |21[ | we find 
the following approximations 


ix?> ^ E 

n=l 

oo 

ixt > ~ E 

n = 2 
oo 

ixt+> ~ E 


(c+ + dpe 




2 V 2 


{c+ + d+)e 




2 V 2 

‘ + (e+-d+)e 


Pn-l|n - 2), 


-Pn-i\n - 1>. 


In order to simplify these expressions we perform a Taylor 
expansion in the frequencies LJn around fi -|- 1 as 


‘^n/g 


\J 4n -|- 2 -f 2 


n — n — 1 
V4fi, -|- 2 


^ (n-n- 1 )'^ 
(4h -h 2)‘3/2 ■ 


(24) 


The previous second order expansion is valid provided the 
third order contribution multiplied by gt is negligible. This 
imposes the restriction on the interaction time 


( 4 h-|- 2 ) 5/2 


nig. 


(25) 


For the rescaled time t = gi/7r\/4fi + 2 used in the main text 
this implies r <IC \/h/27r. In this approximation the field states 
can be written as 


I 0\ -i4,i<^++dp\at -l)-(c+-dp\a+ -\) + 2d^\a) 

2 v *2 

I 1\ (''++‘ij)l‘''rd>-(c+-<ij)|c«+,l>-2d-|c«> 

IXU-e 

('^++d+)]a~ ,0) + (c^-d-^}\a^ ,0) 

IXt / ~ 2 ’ 

with 


(26) 


, ± , vA a" ±i (^+2 

’/) = E ® 

wn\ 


n = 0 

±ij2iTTi ± 


n-n+j 


a*. 


|n> 




(27) 


T 



Fig. 6 Top (bottom) figure: Real part of the overlap (q:|o^} 
((—o|a:)")) as function of the rescaled time t. The red curve 
was evaluated using the exact expression in the first line of 
Eq. ( |28[ ) and the black narrow line corresponds to the ap¬ 
proximation given in Eq. < |29[ ). The mean photon number is 
h = 12.16. 


j S { — 1,0,1} and u = ^Ah -|- 2. Furthermore, the states \cy^) 
are dehned by Eq. ( |10[ l. We neglected the contribution of j in 
the quadratic term of the exponent in Eq. ( |27[ ). This can be 
justified given the fact a Poisson distribution with high mean 
value is almost symmetrically centered around its mean with 
variance equal to its mean. This implies that the maximal 
relevant value in the quadratic term is given by 


max 


(n- n-h j)2 

V 



jj_ 

2n 


which shows that the contribution of j = —1, 0,1 to this term 
is negligible for 1. 

Finally, using the approximations of Eq. ( |26[ ) and \27\ in 
Eq. ( |20[ ) and separating the atomic states accompanying to 
the photonic states |«r ) and la) yields the result of Eq. (HI- 


B Evaluation of (a|a^) and (—a|a^) 


In this appendix we investigate the overlaps between the field 
states laj/) and |iba) defined in Eqs. ( |10[ | and respectively. 
Using the index j S {—1,1} one can write a single expression 
for the four overlaps as 


(ja\Q^) = —-^e 

'.J I T I ^ 

n = 0 

p±i(n+l)27rT °° 




(28) 


yJl'KU 


^ g±» 7 rn( 2 r+V)- “E"’ {n-nf 


n= — 00 
p±i(n+l)27rT °° 


vm 


/ • \ 2 - / \ 2 


n = —OO 


In the second line we have approximated the Poisson distribu¬ 
tion by a normal distribution and we have extended the sum 
to — 00 . These approximations are valid in the limit h ^ 1. In 
the third line we have used the Poisson summation formula 
BZI which in the case of a Gaussian sum can be expressed as 

OO 

^ ^i2nun — s{n — 7i)^ _ 

n= —00 




i 27 vn(n-\-u) — (n + u)^ 


n= — 00 
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with Re[s] > 0. The last expression in Eq. ( |28[ l involves a 
summation of Gaussian terms with variance (l + 7r^T^)/47r^fi. 
This variance is very small provided the condition 4ri 3> is 
fulfilled. If this requirement is met, there exists a dominant 
contribution in the summation that corresponds to the value 
of n where |t + (1—j)/4±n| achieves its minimum value. This 
minimum can be evaluated as 

/j(t) = frac (r + 


where frac(x) denotes the fractional part of x. By considering 
only the dominant term of the last summation in Eq. ( |28[ ) one 
can find the following approximation of the overlap between 
field states 


g±*27r[n/j (T) + (n+l)T] 

~ —- ;■ , . —e~ 

V-1- =t I'KT 




(29) 


with j S { — 1,1}. This result for r = 1/2 and j = —1 has been 
rewritten in polar form in Eq. of the main text, where 
we used that /_i(l/2) = 0. In Eq. ( |13| we have used that 
/i(l/2) = —1/2. In the top panel of Fig. |^have plotted the 
real part of the overlap {o|a^} as a function of the rescaled 
time T. The evaluation of the exact expression is shown in 
red and the approximation in black. The collapse and revival 
phenomena are well described by the approximation of the 
overlap in Eq. | |29[ ). Similar treatment to describe the collapse 
and revival phenomena in the Jaynes-Cummings model has 
been presented in Ref. 1481 . In the bottom figure of Fig. we 
have plotted the real part of the overlap {—o|a^). 


References 

1. H.-J. Briegel, W. Diir, J. I. Cirac and P. Zoller, Phys. 
Rev. Lett. 81, 5932 (1998). 

2. W. Diir, H.-J. Briegel, J. I. Cirac and P. Zoller, Phys. 
Rev. A 59, 169 (1999). 

3. C. H. Bennett, G. Brassard, S. Popescu, B. Schumacher, 

J. A. Smolin, and W. K. Wootters, Phys. Rev. Lett. 76, 
722 (1996). 

4. D. Deutsch, A. Ekert, R. Jozsa, C. Macchiavello, S. 
Popescu, and A. Sanpera, Phys. Rev. Lett. 77, 2818 
(1996). 

5. M. Zukowski, A. Zeilinger, M. A. Horne, and A. K. Ekert, 
Phys. Rev. Lett. 71, 4287 (1993). 

6. N. Sangouard, C. Simon, H. de Riedmatten, and N. Gisin, 
Rev. Mod. Phys. 83, 33 (2011). 

7. P. van Loock, T. D. Ladd, K. Sanaka, F. Yamaguchi, K. 
Nemoto, W. J. Munro, and Y. Yamamoto, Phys. Rev. 
Lett. 96, 240501 (2006). 

8. T. D. Ladd, P. van Loock, K. Nemoto, W. J. Munro, and 
Y. Yamamoto, New J. Phys. 8, 184 (2006). 

9. P. van Loock, N. Liitkenhaus, W. J. Munro, and K. 
Nemoto, Phys. Rev. A 78, 062319 (2008). 

10. D. Gonta and P. van Loock, Phys. Rev. A 88, 052308 
(2013). 

11. J. Z. Bernad and G. Alber, Phys. Rev. A 87, 012311 
(2013). 

12. J. Z. Bernad, H. Frydrych, and G. Alber, J. Phys. B 46, 
235501 (2013). 

13. J.M. Torres, J.Z. Bernad and G. Alber, Phys. Rev. A 90, 
012304 (2014). 

14. C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. 
Peres, and W. K. Wootters, Phys. Rev. Lett. 70, 1895 
(1993). 


15. C. H. Bennett and S. J. Wiesner, Phys. Rev. Lett. 69, 
2881 (1992). 

16. E. Knill, R. Laflamme and G. Milburn,Nature 409, 46 

( 2001 ). 

17. T. B. Pittman, M. J. Fitch, B. C Jacobs, and J. D. Fran- 
son, Phys. Rev. A 68, 032316 (2003). 

18. W. J. Munro, K. Nemoto, T. P. Spiller, S. D. Barrett, P. 
Kok, and R. G. Beausoleil, J. Opt. B: Quantum Semi¬ 
class. Opt. 7, 135 (2005). 

19. Y.-H. Kim, S. P. Kulik and Y. Shih, Phys. Rev. Lett. 86, 
1370 (2001). 

20. C. Schuck, G. Huber, C. Kurtseifer, and H. Weinfurter, 
Phys. Rev. Lett. 96, 190501 (2006). 

21. T. Pellizzari, S. A. Gardiner, J. 1. Cirac and P. Zoller, 
Phys. Rev. Lett. 21, 3788 (1995). 

22. S. Lloyd, M. S. Shahriar, J. H. Shapiro and P. R. Hem- 
mer, Phys. Rev. Lett. 87, 167903 (2001). 

23. F. Schmidt-Kaler, H. Hffner, M. Riebe, S. Guide, G. P. 
T. Lancaster, T. Deuschle, C. Becher, C. F. Roos, J. Es- 
chner, and R. Blatt, Nature 422, 408 (2003). 

24. L. Isenhower, E. Urban, X. L. Zhang, A. T. Gill, T. 
Henage, T. A. Johnson, T. G. Walker, and M. Saffman, 
Phys. Rev. Lett. 104, 010503 (2010). 

25. C. Nolleke, A. Neuzner, A. Reiserer, C. Hahn, G. Rempe, 
and S. Ritter, Phys. Rev. Lett. 110, 140403 (2013). 

26. M. Tavis and F. W. Cummings, Phys. Rev. 170, 279 
(1968). 

27. B. Casabone, A. Stute, K. Friebe, B. Brandstatter, K. 
Schiippert, R. Blatt, and T. E. Northup, Phys. Rev. Lett. 
Ill, 100505, (2013). 

28. B. Casabone, K. Friebe, B. Brandsttter, K. Schppert, R. 
Blatt, and T.E. Northup, Phys. Rev. Lett. 114, 023602 
(2015). 

29. R. Reimann, W. Alt, T. Kampschulte, T. Macha, L. 
Ratschbacher, N. Thau, S. Yoon, and D. Meschede, Phys. 
Rev. Lett. 114, 023601 (2015). 

30. S. Nufimann, M. Hijlkema, B. Weber, F. Rohde, G. 
Rempe, and A. Kuhn, Phys. Rev. Lett. 95 173602 (2005). 

31. R. J. Glauber, Phys. Rev. 131, 2766 (1963). 

32. A. Perelomov, Generalized Coherent States and Their Ap¬ 
plications (Springer-Verlag Berlin Heidelberg 1986). 

33. C. E. A. Jarvis, D. A. Rodrigues, B. L. Gyorffy, T. P. 
Spiller, A. J. Short, and J. F. Annett, New J. Phys. 11, 
103047 (2009). 

34. D. A. Rodrigues, C. E. A. Jarvis, B. L. Gyrffy, T. P. 
Spiller and J. F. Annett, J. Phys.: Condens. Matter 20 
075211 (2008). 

35. M. S. Kim, J. Lee, D. Ahn, P. L. Knight, Phys. Rev. A, 
65 040101(R) (2002). 

36. T. E. Tessier, 1. H. Deutsch, A. Delgado, and 1. Fuentes- 
Guridi Phys. Rev. A, 68 062316 (2003). 

37. J. H. Eberly, N. B. Narozhny and J. J. Sanchez- 
Mondragon, Phys. Rev. Lett. 44 1323 (1980). 

38. W. P. Schleich Quantum Optics in Phase Space (Wiley- 
VCH, Weinheim, 2001). 

39. K. Vogel and H. Risken, Phys. Rev. A 40, 2847 (1989). 

40. J. M. Torres, E. Sadurni, and T. H. Seligman, J. Phys. 
A 43, 192002 (2010). 

41. J. Gea-Banacloche, Phys. Rev. A, 44 5913 (1991). 

42. M. Khudaverdyan, W. Alt, 1. Dotsenko, T. Kampschulte, 

K. Lenhard, A. Rauschenbeutel, S. Reick, K. Schrner, 
A. Widera and D. Meschede, New J. Phys., 10 073023 
(2008). 

43. S. Brakhane, W. Alt, T. Kampschulte, M. Martinez- 
Dorantes, R. Reimann, S. Yoon, A. Widera, and D. 
Meschede Phys. Rev. Lett., 109 173601 (2012). 








Unambiguous atomic Bell measurement assisted by multiphoton states 


11 


44. J. M. Raimond, M. Brune, and S. Haroche, Rev. Mod. 
Phys., 73 565 (2001). 

45. K. M. Birnbaum, A. Boca, R. Miller, A. D. Boozer, T. 
E. Northup, and H. J. Kimble, Nature 436, 87 (2005). 

46. C. Wittmann, M. Takeoka, K. N. Cassemiro, M. Sasaki, 
Gerd. Leuchs, and U. L. Andersen, Phys. Rev. Lett., 101 
210501 (2008). 

47. R.E. Bellman, A Brief introduction to theta functions 
(Holt, Rinehart and Winston, New York, 1961). 

48. A. A. Karatsuba, and E. A. Karatsuba, J. Phys. A, 42 
195304 (2009). 



